ADJUSTMENT COEFFICIENT FOR RISK PROCESSES IN 
SOME DEPENDENT CONTEXTS 
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Abstract. Following [18], we study the adjustment coefficient of ruin 
theory in a context of temporal dependency. We provide a consistent 
0^ ' estimator of this coefficient, and perform some simulations. 
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Adjustment coefficient w for risk processes may describe tlie behavior of 
ruin probability. Several results for sums of i.i.d. claims exist: in [T2], H.U. 
Gerber gave an exact formula for finite time ruin probabilities involving the 
adjustment coefficient w, [19] provide a consistent estimator of w, V. Mam- 
misch [15] gave a necessary and sufficient condition for the existence of w. 
^ ' In dependent contexts, let us cite H.U. Gerber [13] for auto-regressive pro- 

C/^ . cesses, [2] for an extension to ARM A processes and [3l 0] for the study of 

(-H I the adjustment coefficients in Markovian environments . The main objective 

of the parper is to provide a non parametric estimation of the adjustement 
coefficient introduced in [18] in dependent contexts. We give a general de- 
pendent context (weak temporal dependency in the sense of [^) for which 
our estimator is consistent. 
^ , The paper is organized as follows : 

CN I • Section [T] contains the definitions and elementary properties of weak- 

' dependent processes as well as adjustment coefficient. To make 

short, w^, the independent coefficient, will be the adjustment coeffi- 
cient if the process is i.i.d. while w'^ will be the adjustment coefficient 
Q ■ of a dependent sequence. 

0^ i • In Section [21 we prove that w'^ may be seen as a limit (for r oo) of 

' independent coefficients wl- We also provide some general examples 

for which the adjustment coefficient w'^ may be defined. 

• Section [3] is devoted to the estimation of coefficients and w'^ and 
^ • contains the main results : we construct consistent estimators (see 

, Theorems 13.31 13.51 and I3.10p . Note that in |2j , an estimation of w'^ is 

given for ARMA processes which is based on the estimation of the 
ARMA parameters. Our procedure is completely non parametric. 

• In Section m we provide some simulations. 

1. Setting 

We consider {Yn)neN a sequence of random variables and Ru the event 
{Yn > u for some n > 1}. Yn is interpreted as the value of the claim surplus 
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process of a company at the end of the year n. is understood as the ruin 
event for an intial reserve u > 0. We could write 

n 

where Xi is the gain/loss of the company during the year i. 

1.1. Weak dependent processes. This last decade, Doukhan and al. ([71 
El El [H]) have developed new dependence definitions that both extend clas- 
sical probabilistic definitions (as q or <1> mixing) and are satisfied for several 
useful models (like ARM A or ARCH) that are neither a nor $ mixing in 
the standard way. Roughly speaking, in the classical probabilistic definitions 
of mixing, the functions / and g in the definition below (jl.ip belong to the 
whole class of square integrable functions. 

Define (see [7]) for a real valued or vector valued process {Xt)tm-, 
(1.1) eW-supl^""'^'^- 

c{f,g) 

where the supremum is taken over multi-indices i = (ii, . . . , i^), j = (ji, . . . , j^) 
such that: 

k < ■ ■ ■ < iu < iu + k < ji < ■ ■ ■ < jy 

and all functions / : M" — > M,g : W^' — > M are bounded and Lipschitz 
functions, with respect to the distance : 

p 

d{x,y) = ^\xi- yi\, X = {xi, . . . ,xp), y = {yi, . . .,yp). 
1=1 

Remark. We could replace the space of Lipschitz functions by other spaces 
of regular functions (differentiable functions, functions of bounded variation 
...), see [17] for a general condition of convenient functional spaces. 

We define a notion of weak dependence according to the function c{f,g). 

Definition 1. Consider the following functions c{f,g), defined for 

/ : M" — > M and g : — > R bounded and Lipschitz functions, lip(/) 

is the Lipschitz coefficient of the function /. 

(1) c{f,g) = vWfWoo^ipig), we say that the sequence (Xj)jgN is 0-weakly 
dependent if the corresponding mixing coefficients sequence (e(A;))fcgN 
is summable. 

(2) c{f,g) = 'ulip(/)||5||oo + t'||/||oolip(5), we say that the sequence 
(Xj)jgpj is r/-weakly dependent if the corresponding mixing coeffi- 
cients sequence {e{k))k£N is summable. 

This class of dependent processes is very rich and enjoy lots of nice prob- 
abilistic properties. Let us remark that weak dependent processes need not 
to be stationary. 

For completeness, we recall the definitions of a, ^ and ^ mixing. ^ mixing 
may be defined in a formalism close to that of Definition [1] while for a and 
<I> mixing, it is not clear that the same formalism is meaningfull. 
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Definition 2. Consider mixing coefficients e{k) defined by Equation (jl.ip 
wliere tlie supremum is taken over functions 

/ : W — > M and g : W — > M in and c{f,g) = ||/||i||5||i. 
We say tliat the sequence (Xj)jgN is ^ mixing if tlie corresponding mixing 
coefficients e{k) are bounded. 
a and $ mixing coefficients are defined as : 

a{U,V)= sup \¥{UnV) -F{U)F{V)\, 
ueu, vev 



^{l(,V)= sup 

ueu, vev 



nunv) ^j^j 



F{U) 

A process {Xt)^i is a (resp. mixing is the coeficients 

ax (r) = sup a{a{Xt, t<i), a{Xt, t>i + r)), 



resp. ^x{r) =sup$(cr(Xt, t < i), a{Xt, t > i + r)) 



go to 0. 



1.2. Adjustment coefficient. In the classical i.i.d. (i.e. the Xi are i.i.d. 
random variables) model of ruin theory, the adjustment coefficient t/; > is 
defined as the unique positive solution of A(w) = with 

X{t) = logE[exp(tXi)] 

assumed to be well defined. Mammischt ([15]) gave a necessary and sufficient 
existence condition for w. The importance of the adjustment coefficient is 
revealed by the exact formula due to Gerber ([I2]) : let T be the ruin time 
(T = inf{A;GN / Yk > u}), 

^—wu 

¥{Ru) = P(r < oo) 



and the famous de Finetti bound follows : 

We shall focus on the following asymptotic result also due to Lundberg : 

,. logP(^J 

hm = —w. 

As already mentioned above, several attempts to extend these results to 
dependent and/or non stationary settings have been proposed. We wish to 
give a general dependent setting in which such an asymptotic result holds, 
as well as provide a consistent estimator to the adjustment coefficient in this 
dependent context. Our approach does not require a precise knowledge of 
the dependence structure, nor on the law of (Arj)jgN but only an information 
on the speed of mixing (given by Equation (II. ip ). 

Following [18] we assume that : there exists to > such that for all < t < 
to, 

(1.2) c(t) = lim exists. 

Also, there exists < t < Iq such that c{t) = 0. 
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We shall provide a sufficient condition that implies existence and unique- 
ness of a positive solution to c{t) = 0, provided that (|1.2p is satisfied. We 
shall denote w'^ this unique solution. We shall also denote by it;* the unique 
positive solution to X{t) = with 

X{t) = logE[exp(tXi)] . 

Of course, if the sequence Xi is i.i.d. then = w'^. 

1.3. Existence condition. In Mammisch (jl5j) it is proven that exists 
(and is unique as a positive solution to X{t) = 0) if and only if the following 
three conditions are satisfied. We shall denote by (E) these three conditions. 

(1) E(Xi) < 0, 

(2) F{Xi > 0) > and 

(3) either a < oo and E(e°^i) > 1 or a = cxD where 
(1.3) a = sup{t > 0, E(e*-^i) < oo}. 

This condition, together with the weak dependence assumption is sufficient 
to get consistency and asymptotic normality of the estimator of . In order 
to get existence and uniqueness (as a positive solution to c{u) = 0) of w'^, 
we shall need some additional conditions. 

Property 1.1. Assume that the limit c{u) is well defined on [0, no[, in 
particular, for all t < uq, for n large enough, E(e*^") < oo. Assume that (E) 
is satisfied for < a < uq, and 

(1) Cn{t) exists for all n and < t < a, 

(2) for large enough n, ¥(Yn > 0) > 

(3) if a < oo then for n large enough, lim E(e*^") > 1, 

i— »a— 

(4) c'(0+) < or equivalently, 3t > such that c{t) < 

then there exists a unique positive solution to c{u) = 0. This solution is 
denoted by w'^. 

Proof. We adapt Mammish's arguments. Recall that any convex function 
is continuous and admits left and right derivatives on all points where it 
is defined. Moreover, if is a sequence of convex functions defined on 
[0,^0 [ and converging to / on [0,iio[ then / is a convex function and the 
convergence is uniform on any compact subset of [0, uq[. 
The function c is the limit on [0,Uo[ of convex functions 

Cn{t) = ilogE(e*^"). 

n 

Thus c is a convex function with c(0) = and we assume that c'(0+) < 
(which is equivalent to 3 i > such that c{t) < by convexity). 

• If a < oo then, we assume that for large enough n, oo > E(e''^") > 1 
thus oo > c(a) > . Since c is continuous (because it is convex) , we 
deduce that there exists w > such that c{w) = 0. This solution is 
unique because of the convexity of c. 

• If a = +00, then c{t) is well defined for any t G M"*". Because 
¥{Yn > 0) > for n large enough, we have that for large enough n, 
lim E(e*^") = +oo. As a consequence, we have that lim c{t) > 0. 

t— ►oo t— ►oo 



ADJUSTEMENT COEFFICIENT IN SOME DEPENDENT CONTEXTS 



5 



Since we assume that c{t) < for some t, the convexity of c then 
imphes that there exists to such that c'(t±) > for any t > Iq. 
We deduce that c{t) > for t large enough. Since c is continuous 
(because it is convex), we deduce that there exists w > such that 
c{w) = 0. This solution is unique because of the convexity of c. 

□ 

Remark. The condition ¥'{Yn > 0) > is necessary because if it exists u; > 
such that c{w) = then by convexity, either there exists t > such that 
c{t) > or c{t) = +00 for all t > w which implies that E(e*^" ) > 1 for large 
enough n which implies ¥{Yn > 0) > for large enough n. 

2. Limit result and examples 

2.1. Asymptotic behavior for ruin probability. In [TH], it is proven 
that if the adjustment coefficient w"^ exists then it describes the asymptotic 
behavior of the ruin probability in the following sense: 

(2.1) Ita ^ 

U—^00 U 

As a consequence to (|2.ip . we obtain that if it exists, w'^ is the limit of the 
adjustment coefficients of the sequence (yn)neN- 

Corollary 2.1. Assume that the hypotheses of Property are satisfied. 
For large enough n, there exists a unique Wn > such that 

E(e"'"^") = 1 

and w'^ = lim Wn . 

n— >oo 

Proof. The existence of Wn follows from [15j : Yn satisfies the existence hy- 
potheses of Mammish for n large enough. Applying Markov's inequality, we 
get for all K > 0, 

-logP(yn >^^) 

Wn < . 

U 

Then, (|2.1|) implies that any limit point w of the sequence Wn verifies : 
w < w'^. The convergence of functions c„ to c is uniform on [0,w'^], thus 
we have that c{w) = 0, so that either w = or w = w'^. Now, cannot 
be a limit point of the sequence Wn because otherwise we would have that 
c{t) > for all t > which contradicts the hypotheses of Property 11.11 We 
conclude that Wn ■w'^. □ 

Let us give some examples for which the function c{t) is well defined. 
We recall the following result on approximate sub additive sequences due to 
Hammersley |14j . 

Lemma 2.2. Assume h : N — > M be such that for all n,m>l, 

h{n + m) < h(n) + h{m) + A(m + n), 
with A a non decreasing sequence satisfying : 

(2-2) E^^<- 
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h(n) 

Then, A = lim exists and is finite. Moreover, for all m > 1, 

n— >oo n 

Mm) _ A(m)^4 y A(r) 
~ m m ^-^ rir + 1) 

r=2m ^ ' 

Of course, for A(r) = 0(1), then (j2.2p is satisfied. Lemma 12.21 asserts 
that A(r) may go to infinity but not too fast. 

2.2. ^'-mixing processes. According to Definition El we consider the fol- 
lowing classical ^-mixing condition : 

(2.3) >i,w = sup^°-<^'-^- ,;^-,?-?"^--'-^'-"<o°. 

Iilblli 



where the supremum is taken over functions /, g a and over multi-indices 
i = (ii, . . . , i„) and j = (ji, . . . , j,,) with n < • • • < i„ < + A; < ji < • • • < 
jy. The sequence if ^'-mixing is bounded. 

Proposition 2.3. Assume that for t G [0, a[, for all n G N, E(e2*'5") < oo 
and {Xn)n&n is a ^-mixing process then 

lim -logE(e*'^") 

n— >oo n 

exists for any t G [0, a[. 
Proof. Using ()2.3p . we have : 

|lE(e*'^"+'") -E(e*^")E(e*^'")| < ^'(l)E(e*'^")E(e*^'"). 
We conclude the proof by using Lemma 12.21 □ 

Examples of ^'-mixing processes are finite state Markov chains of any 
order but also Variable Length Markov Chains (VLMC) on a finite state 
(see Lemma 3.1 in jllj). 

Even if ^' ($, a)-mixing processes are often used in probability theory, lots 
of useful processes (like ARMA processes) are not ^ (<I>, a)-mixing. In the 
following two subsections, we provide a class of r/ mixing processes for which 
the function c{t) is well defined provided the mixing is sufficiently fast and 
the variables Xi are almost surely bounded. This condition is close to the 
one used in [1] by Brie and Dembo for one other class of mixing processes 
(namely a mixing processes). Then we prove that if the sequence has some 
structure (here, the sequence is a Bernoulli shift) then the condition on the 
speed of mixing may be weakened. We refer to [7] for other examples of 9 
and rj weakly mixing processes (including ARMA and ARCH processes). 

2.3. Super mixing processes. We prove that if the process (Xj)jgp^ is ry 
weakly dependent (recall Definition [T]) with dependence coefficient e(n) = 

O (^6'"(''"")'') with < 6* < 1 and /3 > 1 then the function c{t) is well defined 

provided < Ma.e. A rj weakly dependent process with dependence 

coefficient e{n) = O (e^^^i"")") will be called a super mixmg process. 
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Proposition 2.4. Assume {Xn)neN is a rj weakly dependent process with 
mixing coefficients e(n) = 0(e~'^"''^^°")'^) with c > 0, P > 1. Moreover, 
assume that there exists M > such that \Xi\ < M a.e. Then the sequence 
c{t) is well defined on M. 

k 

Proof. For any < j < k, let Sj = '^Xi, because \Xi\ < M a.e., 

g-Kfc-i)M <jE(gt5^*) <gt(fc-j)M_ 

Also, for any j < i < j, 

Remark that the function : x i— > e*^' is bounded above by e*^'^ and has 
Lipschitz constant te**^ for x € [—M,M]. 

Fix an integer < r < max(n, m). Firstly, assume that n < m and, using 
the definition of r] weak dependence, we get: 

E(e*'^"+'") = E(e*^"e*'^"+i") 

< e''-^ (E(e*^" ) • E(e*^"i"+i ) + e(r) (?i + m)te*'^^e*("^-^)^ ) 

< e2*^*'^E(e*^")E(e*^'") + (n + m)e(r)te2*("+™)*^E(e*^")E(e*^™) 

< E(e*^")E(e*^'")(e2*^*-^ + (n + ?n)e(r)te2*("+'")*^). 

We conclude the proof by choosing r = 0{j^^^^^^^^ ),l < K < (3 and applying 
Lemma 12.21 with a function: 

= O (^) . 

If n > m, the proof is similar but uses the decomposition : 

Q cm. _|_ Q-m+r , Qn+m r— i 

>Jn+m — ^1 -I- '-'m+r+l' '-' 

2.4. Bernoulli shifts. Causal Bernoulli shifts are processes defined as: 

Xn = H{(n-j, j G N) 

with H a measurable function and (^n)nGZ an i.i.d. process. We shall assume 
the following regularity condition on H : define the continuity coefficients 

dn = \\ sup \H{^n-i,i^'^)-H{^n,---,Cl,Uo,U-i,...)\\\oo 
u={uo,U-i,...) 

and assume that the sequence dn is summable. By adapting the arguments 
of [7] we may prove that such a process is 0-dependent with mixing coefficient 
6{n) = dn- 

Proposition 2.5. Assume that {Xn)neN is a Bernoulli shift satisfying the 
summability condition for the continuity coefficients dn- Then the sequence 
lnE(e*'^") satisfies the hypotheses of Lemma \2.^ and thus c(t) is well defined. 
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Proof. We fix a sequence of real numbers u = {uo,U-i, . . .) and we write 

j=i 

e 

j=i 

e 



(2.4) 



'^H{Cj, . . . . . .) - H{Cj, . . . ,^i,uo,n_i, . . .) . 



Using the stationarity of (^n.)neN) we have that 

e e-i 
4 < Yl d3-i = Edj- 

j=i j=0 

Now, C/f and are independent random variables and thus, 

E(e*^"+'") < e*(Er=i*+Er.i<iO]E(e*^r)E(e*^"^"). 
Applying once more ()2.4p . we get: 



E(e*'^"+'") < exp 
If we denote D = dj , we get 



max(n,m) 

4i Yl 



E(e*^")E(e*'^" 



E(e*^"+'") < e^*'^E(e*'^")E(e*'^'") 
and we conclude by applying Lemma l2.2i 



□ 



Remark. In the above proposition, we could replace the hypotheses of summa- 
bility of the sequence (dj)jgN by the summability of 

, \. with A r = > d,-. 

Lemma 12.21 would apply as well and the limit c{t) would be well defined. 

In [7], it is mentionned that stationary ARMA processes are examples of 
Bernouilli shifts. Non linear autoregressive processes may also be examples 
of Bernouilli shifts. Nevertheless, in order to satisfy that the dj are well de- 
fined, it requires that the innovation is bounded. We claim that Proposition 
12.51 remains true for some unbounded Bernouilli shifts but we where unable 
to prove it. 
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3. Estimation 

3.1. Definition of estimators. In this section, we assume that the se- 
quence (Xj)j>i is stationary and that the hypotheses of Proposition 11.11 are 
satisfied. For r G N, the functions K^e''^^) and E(e*^'') may be estimated by 
their empirical moment versions: for A; G N, 

1 

i=l 
i=0 

r 

where = Xj^jj-. Then we define as the positive solution to 

i=i 

logmfc(t) = and Wr as the positive solution to ^ log MJ^{t) = 0. We shall 
prove that w'' is a consistent estimator of and there exists an r = r(k) 
such that Wr = w'^ is a consistent estimator of w'^. We shall also prove that 
they satisfy a central limit theorem. 

Before stating and proving our main results on the asymptotic properties of 
the estimators and Wr, we prove that exists almost surely (the proof 
for Wr will be done later because it requires some weak dependence property 
for (Z[')jgN which is proven in Lemma l3. 40 . 

Proposition 3.1. Assume that the sequence {Xi)i^^ is rj or 6 weakly depen- 
dent and satisfies Condition (E). Then, exists eventually almost surely 
as k — > oo. 

Proof. We begin by noting that the or rj weak dependence implies ergod- 
icity, see for example [6]. Following Section fl.3l we have that exists and 
is unique if and only if, 

(1) iEti^^<o, 

(2) {i = 1, . . . , k / Xi > 0} is not empty. 

Mammisch's third condition is satisfied with a = oo because here, expec- 
tations are finite sums. The two above conditions are eventually almost 
everywhere satisfied because of the ergodic theorem. □ 

3.2. Asymptotic properties of Asymptotic properties of the estima- 
tors w'' and are done by using the same approach as the one used to prove 
results on asymptotic properties for M-estimators in a parametric context 
(see [20] section 5). It is known (see [HllTj) that the process (Aj)jgN satisfies 
a central limit theorem with asymptotic variance 

oo 

r2 = ^Cov(Ao,Ai), 

i=0 

provided that the sequence Cov(Ao, Aj) is summable. 

We obtain a central limit theorem for fhk{t) by proving that the sequence 
^gtx„-j^^j^ is also weakly dependent. Then we prove a central limit theorem 
for the M-estimator 

Let us recall the following results from [7] (Theorem 7.1 and Section 7.5.4). 
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Theorem 3.2. Let (Z„)„gp^ be an r]-weakly dependent sequence with e{n) = 
0(n~2"«^) for K > 0. Then, 

r2 = ^Cov(Zo,Z„), 

n>0 

is well defined and 

As a consequence, we get the following consistency result for as well as 
asymptotic normality. In order to use Theorem 13.21 we first need to prove 
that the sequence {e^^")nef^ is also ?7-weakly dependent with e{n) decreasing 
to zero sufficiently rapidly. 

Theorem 3.3. Assume {Xn)neN is rj-weakly dependent with e{n) = 0{9^), 
< 6 < 1. We have for any t £ [0, uq[ that 

T\t) = ^Cov(e*^«,e*^") 

n>0 

is well defined on [0,^0 [ and 

V^(m„(t)-E(e*^«)) M{0,r\t)). 
converges in probability to w'' and 



with r? 



As already mentioned, we begin by proving that the sequence of random 
variables (e*^")„gN is r/- weakly dependent. 

Lemma 3.4. Assume {Xn)n&n is rj-weakly dependent with mixing coefficient 
e(r). Then, for any t G [0,iio[, the sequence of random variables {e''^'')n£N 
is rj-weakly dependent with mixing coefficient et{r) < 2E(e^*~^'^)"'''^ )e(r) *+« 
with K > such that t + k £ [0, mo[- 

Proof. We follow the proof of Proposition 2.1 in [7]. Let / and g be two 
Lipschitz functions and for Af > fixed, x G M, denote x^*^) = min(x,M). 
Assume (i,j) are multi-indices such that 

ii < ■ ■ ■ <iu<iu+r < ji < ■ ■ ■ < jv, 

and define: 

F{XO = fie'^n , . . . , e'X.^ ) F(^) {X,) = (e<"^ , . . . , e<"\ 
G{X,)=g{e^^n,...,e^x..) G^''^\x,) = (e*^"' , . . . , e^"' ) . 
Then, 

|Cov(F(XO,G(Xj))| < 2||/|UE(|G(Xj)-GW(Xj)|) 

+2||5||ooE(|F(Xi)-FW(Xi)|) 
+|Cov(F(^^)(Xi),G(^-^)(Xj))|. 
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Recall that E(e''-''-') < oo for any s G [0, uq[ and let k > be such that 

t + Ke [0, Mo[, 

then E(|e^*^'^^'''"* I) < oo, using the Markov inequality, we get: 

E(|G(Xj)-GW(Xj)|) < lip(5)^E(|e*^^. -e<"'|) 

k=l 

< 2t>lip((7)e-'^*-^E(e(*+'^)^i). 

Also, since {Xn)n<^N is r/-weakly dependent, 

|Cov(F(^-^)(Xi),GW(Xj))| < (nlip(/) H^IU + ^;hp(5)||/||oo)te*^e(r). 

Finally, we obtain for M > 1, 

|Cov(F(Xi),G(Xj))| < (nlip(/) ||5||oo + ^^hp(g)||/|U) 

x(te*^^e(r) + 2E(e(*+'')^i)e-'^^). 

To conclude, we choose M = — ln(e(r)). □ 

Proof of Theorem \3.3[ Lemma [33] together with Theorem 13 . 2 1 imply that for 
any t G [0,uo[, 

v/^(7n„(0 -E(e*^«)) AA(0,r2(t)) 

provided that the sequence et{r) = e(r)t+« = 0(r~^~") for some a > 0. 
Since we assume that e(r) = 0{9^), < 9 < 1, this condition is satisfied. 
This convergence in law also leads to 

m„(t) "-^ E(e*^«) in probabihty, 

moreover we have that this convergence takes place almost everywhere be- 
cause of the ergodic theorem. 

Now, let us consider the estimator of w^. Following the proof of Lemma 
5.10 in [2^, we have that converges to in probability (this uses the 
convergence in probability of fhk{t)^ the continuity of the map t i— > fhj^[t) 
and the uniqueness of as a positive solution to fhk{t) = 1). The central 
limit theorem follows now from the A method : 



dt 2' ' 9*2 ' 

with w G [min(u;*, u;*), max(i(j*, {(}*)]. Thus, 

^^ + 2^" 

We have that 

= Vfc(E(e'"'^i) -mfcK)) 

and therefore it is asymptotically normal with zero mean and variance 
T'^{w^). Moreover, 
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This quantity converges in probability to E(Xie'" and we have that 

Final 



(w^ — 1/7* ) — goes to zero in probabihty. Finally, we have proven that 

ot'' 



E(Xie"''^i)2 J ■ 

□ 



Remark. We could relax the hypotheses that e{r) decreases to exponen- 
tially fast. It is sufficient that e(r)'+K = 0{r~'^~°') for some a > 0. For ex- 
ample, some intermediate speed of mixing like e{r) = O ^e"^^''^^^''^ , K > 0, 
P > 1 or e{r) = O (6''"°), 0<6'<1, 0<a<lis convenient. 

3.3. Asymptotic properties of w'^. Now, we are interested in the consis- 
tency of w'^. 

Theorem 3.5. Assume that hypotheses of Theorem \3.3\ are satisfied. Then 

there exists a sequence r = r{k) ^ cx) such that w'^ converges in probability 
to w'^. 

Theorem 13 . 5 1 will be proven by rewriting [3]3] for Wj. instead of {D* and then 
by using Corollary 12.11 We only need to prove that the sequence (ZJ)igi^ 
satisfies some weak dependence property. 

Lemma 3.6. Assume that {Xn)n£N is rj-weakly dependent with mixing coef- 
ficient £{k). Then, the sequence is rj-weakly dependent with mixing 
coefficient ez{k) = re{r{k — 1)). 

Proof. Let / and g be two Lipschitz functions. Assume (i, j) are multi-indices 
such that 

ii < ■ ■ ■ <iu <iu + k < ji < ■ ■ ■ < jy. 
Then,Cov{fiZl,...,Zl),g{Z^^^,...,Zl)) 

= Gov ^/(Ajjr+i, . . . , Ajj(r+i), Xjjr+i, • • • , A(j^_,_i)j.) , 

g{Xj-^r+l, • • • , -^ji(r+l)) ^J2r+l> • • • i ^(>+l)r)) 

< reiik - l)r) (^xlip/UfflU + ^||/||oclip/) , 



r 



where ip{xi, ...,Xr,..., Xrk) = I ^ Xj, . . . , ^ 1 . □ 

\j=i 1=1 / 

As a corollary to Lemma [HTUl we deduce that for any r, Wr exists eventually 
almost surely. 

Corollary 3.7. Assume that the hypotheses of Property are satisfied 
and that the sequence {Xi)i^fq is rj-weakly dependent. Then, for any r, Wr 
exists eventually almost surely. 

Proof. This is a direct consequence of Proposition 13.11 and Lemma 13.61 □ 
Proof of Theorem |g.5i From Lemma 13.61 and Theorem 13. 3| we get 
yfcte(i)-E(e*^0) -^M{Q,Tl{t)) 
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with 

r2(t) = ^Cov(e*^S,e*^") 

n>0 

and if we denote by Wr the i-adjustment coefficient of the sequence Zq and 
uir the positive solution to M^{t) = 1, 

Vk{w;.-Wr) -^Af{0,Tl) 

with 

This imphes that 2^ goes to Wr in probabihty, as k goes to infinity. We 
conclude the proof of Theorem 13.51 by using Corollary 12.11 : there exists a 
sequence r{k) '^zl^ qq such that Wr converges to w'^ in probability. □ 

Theorem 13.51 is interesting from a theoretical point of view but it is not 
so useful from a practical point of view. Indeed, it proves that Wr converges 
to w'^ for a sequence r = r{k) but we have no information on how to choose 

r with respect to k. Moreover, provided that — ^ ^^^z^ 2 converging, we 

could obtain a central limit theorem for ^/k{vur — w'^) with limit variance 

= lim — ; "tfttt- This expression of the asymptotic variance in the 

central limit theorem is not useful from a practical point of view. We might 
nevertheless use moment and Bienaime-Tchebitchev inequalities in order to 
get a useful relationship between r and k. 
We apply the following inequality on the order 2 moment. 

Proposition 3.8. Let (Wj)jgN be a sequence of centered random variables 
and 

C,-2= sup Cov{Wt„Wt2). 

t2-tl=j 

n-1 
j=0 



where 5„, = Wi 



1=1 

As a consequence, if Wi is rj-dependent and stationary, with mixing coeffi- 
cient e{k) < C0^, we have that 

E{Sl) < WnMic-^ ^ 

1 — m 

where Mm = E(|Wip). 

Proof. The first part of the proof of Proposition 13.81 may be found in 
(see Lemma 4.6 p. 79). For the second part, we proceed as in the proof of 
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Lemma E3] (see also [TO]) : let wj;^^ = max(mm(Wt, M), -M), (so that 
^ provided that \Wt\ < M). Then, 

CoY{Wt,Wt+j) < Cov(W^f'\vr,(^f) + Cov(W^f'\(t^i+,--T^(|^^ 

< 2Mh{j) + 2M\\Wt - + \\{Wt - ^Vlli 

< 2M'^e{j) + GM^-'^Mm, 
where the last line is obtained by noting that 

Wt-W^^'''\\i = j l\w,\>M\Wt-Wi^'\dF 

< 2 j l^w4>M\Wt\dF 

< 2M^~™Mm. 

Similarly, we obtain 

We conclude by choosing M = ( — ^ | . □ 

ve(j)y 

Proposition 3.9. Assume that (Xj)jgN is an r] weakly- dependent process, 
with mixing coefficient e{k) < C9^, C>0, O<0<1. Then, for any 
t G [0, uq[, such that 3t < uq, for any v > 0, we have: 

4(C7r + 3)E(e3*^'-)2/^ 



M[(t) -E(e*^'-)l > < 



Proof. We apply the Bienaime-Tchebitchev inequality and Proposition [37 
with m = 3. Let 

Following the lines of the proof of Proposition 13.81 we have that the covari 
ance coefficients associated to (Wf(t))£gN are: 

C2,j < M^6i'^{2Cr + 6), 
with Mm = E{e*'^^-). We have used the fact that 

CoY{wi^'\t),wi^'^{t)) < 2Mlre{r{j - 1) + 1) < 2M^rC0i 
We choose m = 3 and apply Proposition 13.81 so that 



E{[kMl{t) - A;E(e*^'-)]^) < 2A;(2Cr + 6)[E(e'^*^'-)]i 



1-9^ 



We conclude by using the Bienaime-Tchebitchev inequality. □ 

This proposition shows that in order to get the consistency for Wr we 
should choose r{k) = o(lnA;). 

Theorem 3.10. Assume that (Xj)jgN is an rj weakly- dependent process, 
with mixing coefficient e{k) = 0{9^), < 6* < 1 and that 3w'^ < uq. Then, 
for r = r{k) = o(lnA;), Wr goes to w'^ in probability. 
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Proof. We have that 
E(e'"-^r) -M[(u;,,) = Ml{wr) -Ml{wr) 



{Wr - Wr) H / dw, 



dt ' 2 J dt"^ 

Tr 

with Ir = [mm{iVr,'Wr),nia'K{'Wr,Wr)], so that 

dMliwr) 1 f d'^Mliw) 



\Wr — 'Wr\ 



|E(e-^^r)-M[K)| 



dt 



+ 2 



dw 



For any L > 0, 



dMl{wr) 
dt 



remark that —^^ > 0. 



dt"^ 



¥{\Wr — Wr\ > u) 



dMliWr) 



dt 



< P(|M[K) - E(e"''-^i )l >uL) + ¥ 



dMliwr 



dt 



<L\. 



Denote a^(t) = E (^^^^ 



|E(e*^'-) and take L = 



Note that ar{wr) > because of the convexity of E(e*^i ) and the fact that 
Wr is the unique positive solution to E(e*^i) = 1. Hence, 



dMliwr 



dt 



< L 



< 



< 



dMliwr) 



dt 



ar{Wr) 



> 



ar{Wr) 



16(C7r + 3)E(y^3g3«;.y.) 



kar{Wry{l-e6) 

(we proceed as in the proof of Proposition 13. 9p . Finally, 

(3.1) ¥{\wr-wr\ >v)< 16(Cr + 3) (]E(e3»..y.)§ i_^E(y,=^e^""-^-)§). 

Remark that ar{wr) = rE(e^'"^'')Cr{wr)- Consider an interval [ui,U2] C 
[0,^0 [ such that w"^ G [ni,ti2], 3u2 < uq, c is non decreasing on [ni,U2], 
c{ui) < 0, c{u2) > (this exists because of the convexity of the function 
c). Since Cr converges uniformly to c and Wr converges to w'^, for r large 
enough, Wr G [ui,U2], and 

I _ gr(c(tii)-e) 

We also have that for r large enough, ar(wr) > . By taking 

U2 - Ui 

r{k) = o(lnA;), we have that Wr goes to w'^ in probability. 



□ 
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4. Simulations 

We conclude by giving some simulation results. We present some models 
for which the adjustment coefficient is computable - namely MA and AR 
linear processes with an innovation following an exponential law. We provide 
also a non linear example. We refer to [5] for non linear and computable 
examples. These examples are also more realistic from an actuarial point of 
view. 

Recall that if follows an exponential law with parameter 6 > then for 

o<t<e, 




Let £i = — c, with c9 > 1. Then, the independent adjustment coefficient 
u)* is the positive solution to : 



The simulation results are summarized below. The graphs represent the 
estimator log mk{t) and -logM^(t) of X{t) and Cr{t) respectively in grey 
and in black. 

4.1. Independent case. We have simulated an independent sample of = 
— c of length 10000 and 9 = 1.2, c = 1. We have computed and w"^: 



r = 6 


w' = 10^^ = 0.38 


w' = 0.36 


w'^ = 0.37 




\ \ \ \ \ I 

0.30 0.32 0.34 0.36 0.38 0.40 



4.2. AR{1) model. We consider the following ^i?(l) model: Xn = aXn-i + 
En- Following [13], we have that w'^ = {l — a)w^. We have simulated a sample 
of length 10000 for 6 = 1.2, c = 1, a = 0.3. Then, 



r = 6 


= 0.45 


w"^ = 0.26 


w'^ = 0.27 
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4.3. MA{1) model. We consider the following MA{1) model: X„ = e„ + 
aSn-i, with 9 = 1.2, c = 1, a = 0.2. Then, w"^ is the positive solution to: 

-tc(l + a) + ln6' - ln(6l - t{l + a)) = 0. 

We have simulated a sample of length 10000, 



r = 6 


= 0.47 


w'^ = 0.31 


w"^ = 0.32 




4.4. A non linear AR{1) model. We consider the following non linear 
AR(1) model (which may be seen as a particular case of Bernoulli shifts, 
see [7]): Xn = aX^_^ + 0.7e„. We have simulated a sample of length 10000, 
with 9 = 1.2, c = 1, a = -0.2. 



6 



0.8 



w 



1.21 



18 H. COSSETTE, E. MARCEAU, AND V. MAUME-DESCHAMPS 



4.5. How to choose the r-parameter ? When performing the estima- 
tion of the w'^ coefficient, we are faced with the choice of the parameter 
r. Fohowing Theorem 13.101 we should take r = o(ln(A:)) but the practical 
choice of r for n given is not clear. We have performed several simulations 
for the independent, MA[1)^ AR[1)^ non linear AR{1) models, for several 
values of r. These experiments tend to show that when r increases, the 
estimator w'^ behaves monotonically in the beginning and then has a more 
chaotic behavior. We propose to choose r as the largest integer for which 
w'^ is monotonic on [0,r]. This is illustrated in the graphs below for several 
models. 



4.5.1. Independent case. We have simulated an independent sample of £i = 
— c of length 10000 and 9 = 1.2, c = 1. Below is represented w*^ for 
r = 1,... ,35, = u;"^ = 0.38. 
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4.5.2. Linear MA{1). We consider the following MA{1) model: Xn = + 
ae„_i, with 9 = 1.2, c = 1, a = 0.3. We have simulated a sample of size 
10000 and represented below u;'^ for r = 1, . . . , 40, w'^ = 0.26. 




4.5.3. Linear AR(1). We consider the following yli?(l) model: X„ = aX„_i+ 
En- We have simulated a sample of length 10000 for 9 = 1.2, c = 1, a = 0.4, 
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4.6. Non linear AR{1). We have represented below w'^ for r = 1, ... ,35 
for the non hnear AR{1) model of section 231 




5 10 15 20 25 30 

1:35 



4.7. On the empirical distribution of w'^. We conclude this simulation 
section with a short study of the empirical distribution of w'^. We have 
performed 100 simulations of a sample of size 10000 in the MA(1) model 
(Section 14. 3p . The mean value of w'^ is 0.317, with standard deviation 0.04. 
The computed value of w'^ is 0.314. The histogram and a Shapiro test 
indicate that the distribution of w'^ is probably asymptotically normal. 
Shapiro-Wilk normality test 
W = 0.9871, p-value = 0.4462 
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